clear;clc
c=12;
k=5.6;
t=0:0.5:600;
F=rand(1,length(t));
v=c*(-log(1-F)).*(1/k);
plot(t,v)
xlabel('t/时间');
ylabel('v/风速');